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Abstract 

Three-magnon problem for exactly rung-dimerized spin ladder is brought up 
separately at all total spin sectors. At first a special duality transformation of the 
Schrodinger equation is found within general outlook. Then the problem is treated 
within Coordinate Bethe Ansatze. A straightforward approach is developed to ob- 
tain pure scattering states. At values S = and S = 3 of total spin the Schrodinger 
equation has the form inherent in the XX Z chain. For S = 1, 2 solvability holds 
only in five previously found completely integrable cases. Nevertheless a partial 
5 = 1 Bethe solution always exists even for general non integrable model. Pure 
scattering states for all total spin sectors are presented explicitly. 



1 Introduction 

Among other gapped ID systems spin ladders were intensively studied during the last 
15 years experimentally, numerically and theoretically (see Refs. in [1]- [3]). The interest 
is accounted by their possibly relation to high temperature superconductivity, variety of 
static and dynamical properties and even an existence of several reliable compounds. 
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In the pioneering paper [3] a spin ladder was suggested as a double spin chain with 
Heisenberg interactions both across and along the chains direction namely rung and leg 
exchanges related to the couplings J T and Ji. It was also pointed that the case, 



has a principle interest because it belongs to the so called rung-dimerized phase in which 
almost all spins are coupled into rung-singlets (rung-dimers). In the purely Heisenberg 
reference model this phase becomes exact only for J/ = 0. However it is always assumed 
that under the condition (1) the physical picture does not change in common. 

Soon it became clear that spin ladder Hamiltonian also admit a term related to di- 
agonal Heisenberg coupling as well as four spin terms [5]. At a first sight these new 
interactions seemed to be complications for a theoretical analysis. However even in [6] it 
was noted that a special linear condition (the Eq. (22) of the present paper) on the former 
and new coupling constants guarantees (for rather big J r ) exactness of the rung-dimerized 
ground state. Besides in this case all one- and two-magnon states also may be obtained 
in explicit form [6], [TJ. 

Unfortunately the rung-dimerization condition (22) has no reliable atomic level inter- 
pretation, so there is no physical reason to postulate it. Nevertheless it seems reasonable 
to suppose that for strong rung exchange any deviations from the exact rung-dimerized 
picture should be small and may be evaluated perturbatively. (In more detail this question 
will be studied in a forthcoming paper.) Under this point of view exactly rung-dimerized 
spin ladders are the best reference models for treating the whole rung-dimerized phase. 

Some static and dynamic zero-temperature properties of exactly rung-dimerized spin 
ladders were studied in a series of papers [7J -[E]- Due to the gap it succeeded to describe 
Raman scattering [7j , magnetic phase transition [8] and (for asymmetric ladders) magnon 
decay [9], [10] utilizing only one- and two-magnon spectrums. Three- magnon problem is 
less actual for the T = physics (see however the papers [H],[T2] devoted to the 5=1 
Haldane chain and 0(3) nonlinear cr-model). 

Advancement into the T > region needs a knowledge of the whole spectrum [3].|16j. 
However such level of clearness may be achieved only for a rather limited list of the so 
called integrable models [3],[l3]-[I7]. The latter besides are significant in heat transport 
phenomena [18J. 

But how to find an integrable model? How it may be identified from a overwhelming 
majority of nonintegrable ones? The most direct way is to express a treating Hamiltonian 
density as a derivative of the corresponding .R-matrix which satisfy the Yang-Baxter equa- 



J r 3> Ju 
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tion. Solvability of this problem is governed by the Reshetikhin condition |17j.[19].|20j. 
If the latter is satisfied for a given local Hamiltonian density then the corresponding R- 
matrix rather exist and may be obtained by an analysis of power series [20] , [21] , [22] or 
by some Yang-Baxterization ansatze [23]. In the present paper we suggest an alterna- 
tive approach based on solvability of the three magnon problem in a framework of the 
Coordinate Bethe Ansatze (CBA) [24"] . 

The essence of the CBA method [T3] , [14j is an assumption that any many-particle 
wave functions is in fact a linear combination of terms produced by multiplications of 
one-particle exponents. Namely, for a rung-dimerized spin ladder the one-magnon wave 
function ip(n) = e lkn [6] is parameterized by a real < k < 2it (the wave number) and 
depends on an integer n (position of the triplet rung). A two magnon wave function 
ip(m,n) (m < n) is linear combination of two exponents e l ( fcim+fc2T1 ) and e^ fe27n+fcin ) [7] and 
so depends on a pair of non equal parameters k\ and k 2 . For a scattering state they both 
are real and one may put 

< k x < k 2 < 2tt, (2) 
while for a bound state they are complex conjugate 

k 2 = h- (3) 

In this light it is seems reasonable to search for representation of multi-magnon wave 
functions as sums of the Bethe exponents. However even a subsequent development of 
this approach to the three-magnon sector dashes on the problem of non integrability. 

In order to reveal an origin of this obstacle let us at first turn back to a two-magnon 
state. Total quasimomentum (wave number) and energy of the latter are the sums 

k = kx + k 2 , E(k 1 , k 2 ) = E magn {kx) + E magn (k 2 ), (4) 

where E magn {k) a single magnon energy. It is significant that under the conditions (2) or 
(3) the mapping 

h, k 2 — > k, E, (5) 

given by (4) is uniqually (up to an exchange k\ ^ k 2 ) reversible. However for three 
magnons the situation is drastically different. Indeed a system of relations 

k = kx + k 2 + k 3 , E = E magn (ki) + E magn (k 2 ) + E magn (k 3 ), (6) 

can define an infinite number of triples (ki, k 2 , k 3 ). As a result a three-magnon wave 
function related to the pair (k, E) should contain in general an infinite number of expo- 
nential terms related to different solutions of the system (6). Evidently such three magnon 
problem is practically unsolvable. 
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The above obstacle may be overcame by existence a first integral (a translationary 
invariant operator commuting with the Hamiltonian) which produces the third condition 
additional to (6). An integrable system has an infinite number of such commuting in 
pairs first integrals and may be solved in all multi-particle sectors. It is significant that 
within the CBA a difference between integrability and non-integrability manifests just at 
the three particle level. As a consequence of this fact one may consider solvability the of 
three-particle problem as an alternative integrability test. 

In the present paper we study three magnon sector of a rung-dimerized symmetric spin 
ladder. At first we briefly analyze the problem in general outlook and only afterwards turn 
to CBA. Motivation of such approach is the following argumentation. Usually CBA is 
treated as a successful ad hoc conjecture which allows to obtain in a rather straightforward 
manner all multi-particle states for a given an quantum integrable model. However the 
reference one is not integrable at general values of coupling constants. As a result (it will 
be shown below in detail) the CBA approach is applicable only in five special integrable 
cases. 

The calculations are performed separately in the sectors S = 0, 1, 2 (the S — 3 sector 
is similar to the S = one) of total spin. At S = (S = 3) the system of equations 
on Bethe amplitudes has a well known form inherent in the XX Z spin chain and so is 
completely solvable for all values of coupling constants. For 5* = 1 and 5* = 2 a complete 
solvability takes place only in the five integrable cases obtained earlier [21] within the 
Yang-Baxter framework. However even in the general nonintegrable case there is a special 
(very complicated) solution in the 5 = 1 sector. Its interpretation remains unclear. 

The plan of the paper is the following. In Sect. 2 we represent the spin ladder 
Hamiltonian in the most tractable form for which the rung-dimerized condition is evident. 
In Sect. 3 we show that the Bethe form of the two-magnon wave function readily follows 
from a straightforward treatment of the Shrodinger equation. In Sect. 4 treating within 
general framework the S = (S = 3) sector we reveal a duality transformation of wave 
function (generalized in Sect. 5 and 6 for S = 1,2) and show that the Bethe Ansatze 
readily follows from the factorized (Fourier) substitution. We also obtain a classification 
(generalized in Sect. 5 and 6 for S = 1,2) of Bethe three-magnon states related to 
complex wave numbers. Pure scattering states obtained within a straightforward approach 
developed in Sect. 4,5,6 are presented in the Appendix. In Sect. 7 we show that the 
revealed CBA solvability is in one to one correspondence with the integrability revealed 
earlier within the Yang-Baxter framework [21]. We also present the corresponding R- 
matrices. In Sect. 8 within CBA we describe action of the S3 permutation group in all 
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total spin sectors. This symmetry as well as duality described in Sect. 5 and 6 is used in 
the Appendix for more compact representation of Bethe states. 

Since the ground state of the model has a simple factorized form we treat it only 
in the infinite volume limit. Analogous approach to the ferromagnetic XX Z chain was 
developed in [T3] . 



2 The spin ladder Hamiltonian 

Before presenting the spin ladder Hamiltonian let us introduce the following local opera- 
tors 

*n = ^(Sl,n - S 2 ,n) ~ i[Sl,n X S 2 ,n], 

*n = ^(Sl,„ - S 2 , n ) + i[S ltn X S 2 ,n], (7) 

(we use the notation ty n instead of more convenient or V&t only in order to avoid 
such rather cumbersome notations as Here Si >n and S 2i „ are local spin operators 

associated with ra-th rung. They may be expressed from and *fr n as follows 



o 2 ,„ - 2 

The representation (8) is similar to the one suggested in [25J but in fact is not identical 
to it. Really the analogs of and \Ev treated in [25] act in an extended vector space. 
That is why for example the "inverse" representation (7) fails for them. 
It may be readily proved that 

Qn\ ^ ni Y^ni Qn\ (9) 

where 

Qn = 2^n' S ra = Sl,n + &2,n, (10) 

Let |0)„ and |1)„ be correspondingly singlet and triplet states associated with n-th 
rung. From (10) follows that 

Qn\0)n = 0, Qn|l)n=|l)n, (H) 
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so the local operator Q n is projector on the n-th rung triplet sector. Then according to 
Eq. (9) the two triples and may be treated as rung-triplet creation-annihilation 
operators. Namely the tripe (a = x,y,z) for which 

*S|0>„ = |1>S, »^|l>i = 0, n\0)n = 0, *l\l) h n = 5 ab \0) n , (12) 

gives the following representation of the total rung-spin 

S a n \l) b n = ie abc \l) c n , (13) 

(e a bc is the Levi-Chivita tensor). Parallel with (12) we shall use the triple 

|1)£=$£|0) B , S*|l>i = j|l>^, j = -1,0,1. (14) 

related to operators 

*f = i=(^±^), = (15) 

It seems reasonable to represent the Hamiltonian density H n ^ n+ i for general spin ladder 
Hamiltonian 

H = J2 H ^n+l, (16) 

n 

in the following form 

MQn + Qn+l) + J 2 (*n " *n+l + *n " *„+l) 
+ JzQnQn+l + J&n ' S n+ i + J$(S n • S n+ i) 2 

+ J 6 (f n + (17) 

Up to a constant this representation is equivalent to the standard one p]-[6] 

H n ,n+i = JrH r nn+1 + J[H l nn+1 + JdHn n+1 + J rr H" n+l + JuH l ^ n+1 + Jdd,H^ n+1 , (18) 
where 

H^ n+1 = — (Sl )7l ■ S2, n + Si jn+ i ■ S2,n+l)? -^n,n+l = &l,n ' &l,n+l + S2,n ' S 2)n +1, 
rd 



Hn,n+1 ~ ^l,n ' S2,n+1 + §2^ ' Sj^+J, H^ n+1 — (Si jn • S2, n )(Sl in +l ■ S2 jn +l); 
jll 

l n,n+l 



n+l — (^l,n • Si ) „ +1 )(S 2) n • S 2 , n +l), n+l ~ (^l,n ' S 2; n+l)(S2,n ' Sl j7l+ i 



and 



Ji — - ^2 J T — 3J rr — Ja — Jddj , 



J 2 = ^(4(Ji-Jd) + Jii-J, 

1/3 Jrri 

J 4 = g( 4 ( J i + J d) + Jll + Jdd), 



J5 — ^ ^JfJ + Jdd 

■h = \[^{Ji- Jd)- J11 + Jdd). (20) 
It was suggested in [5] that only the case 

Jrr = J 11 = — Jdd, (21) 

(or equivalently J5 = 0, Jq = Ji — ^3/2) has a reliable interest. However since spin ladders 
with failed condition (21) also are currently studied [3] and we shall not require it. 
From (9) and (17) directly follows that for 

J 6 = ^ J n - J dd = 4(J, - J d ), (22) 

(triplet-rungs pair creation-annihilation processes are suppressed) there holds 

[H,Q}=0. (23) 

Here the global operator 

n 

according to (11) may be treated as a number operator for rung-triplets. The commutation 
relation (23) results in splitting of the Hilbert space on an infinite sum of eigenspaces 
related to different eigenvalues of Q. In particularly for rather strong J\ the (zero energy) 
ground state of the model has a simple tensor-product form [6] 

|0> = n®l°>«. (25) 

n 

At the same time the physical Hilbert space is subdivided into a direct sum of magnon 
sectors 

00 

n = Y^H m , Q\ Hm =m. (26) 

m=0 

Only this special case (Eq. (22) and rather strong Ji) will be studied in the present 
paper. Additionally we shall imply that J 2 7^ 0. The completely diagonal frustrated 
model related to J2 = or equivalently J d = Ji (in this case the Hamiltonian density (17) 
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may be expressed only in terms of Q n and S„) was studied in details in [26]. Besides one 
may assume that 

J 2 > & Ji > J d . (27) 

Indeed the case J 2 < can be reduced to (27) by use of the following exchange of the 
coupling constants 

Ji <-> J d , Ju <-> J dd , (28) 
related to permutation of spins on all even (odd) rungs. 

3 One- and two-magnon states 

Taking into account (17), (11), (12) and (14) one gets the local formulas 

H n , n+1 . . . |l) n |0)n+l • • • = Jl ■ ■ ■ \l)n\0)n+l ■ ■ ■ + h ■ ■ ■ |0>„|l)„+i ■ ■ ■ , 
#n-l,n...|0> n -l|l>„... = J 1 ...|0) n _ 1 |l) n ... + J 2 ...\l) n ^\0) n ..., (29) 



and 



#n jn+ ie a6c ...|l)^|l)£ +1 ... = e l e abc ...\l) b n \l) c n+1 ..., (a,b,c = x,y,z) 



Here 



and 



Hn,n+1 ■ ■ ■ |l)n|l)n+l • • ■ ~ £ ■ ■ ■ I l/n I Vn+1 • • • > 

c 

n+1 

#„, n+1 ...|l)+|l>+ + i ...= e 2 ...|l)+|l)+ +1 ... (30) 

e s = 2(J X + J 2 A S ), (31) 
J 3 - 2 J 4 + 4J 5 A{J d - 2Ji) + 2J„. + 3 J u 



A 



2J 2 4(J, - J,) 

J3 — J4 + J5 4( J rr — Ji) + J;/ 



2J 2 8(J - J 



d 



J3 + J4 + J5 = 4(2 J d - J t ) + 4J rr + 3 J u 
2 2J 2 8(J,-J d ) ' 1 J 

From (30) may be readily obtained a useful formula 

H n ,n+i . . . \i)M + t ■■■ = (2J1 + j 3 + j 5 ) . . . \mi)Ui ■ ■ ■ + J 4 • • • ii)nii):+i • • • 

+ 5ab(j 5 - Ji) ...\l) c n \l) c n+1 ..., (a,b,c = x,y,z). (33) 



Turning to excitation states we notice that an explicit form of a one-magnon state 



n-1 



\l,k} = J2e lkn ( II ®l°)m)®|l)n®( II (34) 
n m=— oc rrt=n+l 

directly follows from (23) and translation symmetry 

P|l,fc) = e - <fc |l,fc). (35) 

Here P is the translation operator 

P II ®lxW)n = II ®lxW>n+l, XH = 0, 1. (36) 

The corresponding dispersion 

E m agn(k) = 2(Ji + J 2 COS fc) , (37) 

readily follows from (29). 

Since the Q = 2 sector is subdivided on the total spin S = 0,1,2 subsectors we denote 
at once a two-magnon state with total spin S and wave vector A; as |2, S, k). The following 
general representations for the two-magnon states 

|2,0,fc) = ^e^ m+ ")/ 2 a (A;,n-m)...|l)^...|l^..., 

m<n 

\2,l,k) a = e abc J2^ Hm+n)/2 ^(k,n-m)...\l) b m ...\iy n ..., 

m<n 

\2,2,k)+ 2 = J2 eiHm+n)/2a ^ n - m )---\ 1 )t---\^---, (38) 



agree with the rotational and translation (35) symmetries. From Eq. (38) by ". . ." we 
denote an appropriate tensor product of rung-singlets (similar to products in (34)). For 
simplicity the S = 2 sector is represented in (38) by the S z = +2 states. Besides we 
suggest that the reduced wave function as(k,n) should be bounded 

sup as(k, n) < oo. (39) 

n 

The Schrodinger equation for as(k,n) has different forms at n > 1 and n = 1. In the 
former case Eqs. (29) and (30) give 

k 

Uia s (k, n) + 2J 2 cos -[a s (k, n - 1) + a s (k, n + l)] = Ea s (k, n), (40) 
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while in the latter 

(2Jx + E S )a s (k, 1) + 2J 2 cos ^a s (k, 2) = E{k)a s {k, 1). (41) 

It is convenient to rewrite Eq. (41) in the form of Eq. (40) [13] by continuing as(k, n) into 
unphysical region n = 0. Comparing (40) and (41) one conclude that this trick entails a 
Bethe condition 

A s a s (k,l) =cos-a 5 (M). (42) 

The system (40) (considered now for n > 1) and (42) together with the restriction 
(39) allows us to obtain entire as{k,n) in a straightforward manner. Indeed representing 
Eq. (40) in an equivalent matrix form 

a s (k,n + l) \ ( a s (k,n) 

= I , (43) 

a s (k,n) J \ a s (k,n- 1) 
where 

(2k -1 \ E - 4Ji 

10/ 4J 2 cosfc/2 

and taking a5(/c, 1) : ag(fc, 0) from (42) one consequently obtains (up to a constant factor) 
using (43) the rest of as(k, n) at n = 2, 3, ... In following we shall study this problem in 
detail considering separately three regions |k| < 1, |k| > 1 and \k\ = 1. 
For \k\ ^ 1 the matrix J-{k) has two different eigenvalues 



A±(k) = k± Vk 2 - 1, (45) 

related to eigenvectors 

/ A±(k) \ 

e±(«) = I ±V J . (46) 

At \k\ < 1 it is more convenient to use the following representation 

A±(K)=e ±iQ , n = cosq, < q < vr. (47) 
According to (43) a decomposition 



a s (/c, 1) 
a s (k, 0) 



c + e+(K) + c_e-(K), (48) 
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(c± are some coefficients) results in 

= A" (k)c+£ + (k) + A»(K)c_e_(«). (49) 



as(k,n + 1) 
a s (k,n) 



or equivalently 

A; 

a s s catt (k,q,n) = cos- sin - A 5 sing(n - 1). (50) 

The expression (50) obviously agrees with (39) and (as it readily follows from (29) and 
(30)) corresponds to dispersion 

E scatt (k,q) = 4^Ji + J 2 cosgcos|j. (51) 
According to the following formulas 

e ik(m+n)/2 a scatt( k ^ n _ m) = ^ C^^e*^ 1 "^* 2 ^ - C 5i21 e i(/c2m+fcin) ] , 

E sca tt(k,q) = E magn {ki) + E magn (k2), (52) 

where 

- - g = fci < fc 2 = - + g, g = — - — , (53) 

and 

^ = cos^±^-A s e^^)/ 2 , (54) 

one may associate (50) with a scattering wave function of two magnons with wave vectors 
k\ and ki reduced to the center mass frame. 

For |/c| > 1 both the eigenvalues (45) are real. More specifically at ±/c > 1 there 
should be |A T (/t)| < 1 < |A±(/c)| and the representation (49) agree with (39) only for 
c± = 0. According to (42), (46) and (48) in both the cases the eigenvalue remaining in 
(49) is (cos k/2)/A s . So one gets 



E bound (S,k) = 2(2J 1 + J 2 A s + ^-cos 2 ^. (55) 



This solution exists only for 



|A 5 | <cos^ < |A 5 |, (56) 
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and according to the formulas 

^ik(m+n)/2^lK>undff^ ^ m ^ _ gi(fcim+fc2n) 

Ebound(S, k) = E magn (ki) + E magn {k2) , (57) 

where 

h = --iv, k 2 = - + iv, v = \n( — ), (58) 

2 2 \cosk/2/ 

may be associated with a two-magnon bound state wave function reduced to the center 

mass frame. Both (52) and (57) reproduce the Bethe Ansatze calculation presented in [7J. 

For k 2 = 1 the matrix T{k) has an eigenvector £o and an adjoint vector £o 

(59) 

y i i y o i 

for which 

F(k)€o(k) = F(k)£o(k) = «lo(«) + ^o(«)- (60) 

Taking into account (60) and (42) one gets the following decomposition 

as(k, 1) \ k 

= A s Uk) + (cos--kA s )Uk). (61) 
a s (k,0) J 2 

The resulting wave function 

k 

a s (k, k, n) = nK n_1 (cos - - kA s ) + K n A s . (62) 

agrees with (39) only on the appropriate bound of the interval (56), namely for 

k 

cos — = k,As- (63) 

The solution (62) may be obtained from both (50) or (55) in the limit \k\ — > 1. Indeed 
despite the wave function (50) turns to zero at q = 0,ir the ratio a scatt (k,q,n)/ sin q 
remains finite and gives (62) as a limit value. Analogously using the formula (1 + e) n = 
1 + ne + o(e) one can obtain (62) from (55). 



4 S = and S = 3 three-magnon sectors 

Representing at once a S = state in general translationary covariant form 
13,0.*) = e abc e lk(m+n+p ^ 3 b (k,n-m,p-n) 

m<n<p 

...\l) a m ...\l) b n ...\l) c p ..., (64) 
12 



one readily obtains from (29) and (33) the Schrodinger equation for the reduced wave 
function bo(k,m,n) both for to, n > 1 

QJ 1 b (k, to, n) + J 2 [e~ ik/3 b (k, m+l,n) + e ik/3 b (k, m — l,n) 
+e~ ik/3 b (k, to - 1, n + 1) + e ik/3 b (k, to + 1, n - 1) + e~ ik/3 b (k, m,n-l) 

+e ik / 3 b (k, m,n + 1)] = £6 (A;, m, n), (65) 

and m,n — l 

(4Ji + ei)6o(A;, 1, n) + J 2 [ e - i/c / 3 6 (A;, 1, n - 1) + e lfc / 3 6 (A;, 1, n + 1) 

+e ifc / 3 6 (A;, 2, n - 1) + e- ik / 3 b (k, 2, n)] = S6 (A:, 1, n), 
(4Ji + ei)6o(ife, m, 1) + J 2 [e ife/3 & (^, m - 1, 1) + e~ ik / 3 b (k, m + 1,1) 

+e~ ik /%(k, m - 1, 2) + e ik / 3 b (k, m, 2)] = £6 (A;, m, 1). (66) 

Reduction of (66) into (65) results in a system of Bethe conditions, 

2Ai&o(M,n) = e ik/ \(k,0,n) + e- ik/3 b (k,0,n+l), 
2Ai&o(A;, m,l) = e- ife/3 6 (A;,m,0) + e ifc/3 6 (A;,m + 1,0). (67) 

The pair (65) (considered for m, n > 0) and (67) represents the Schrodinger equation for 
bo(k,m,n). It is invariant under the following duality transformation 

V(bo(k, m, n)) = t>o(k, n, to). (68) 

Autodual and anti-autodual solutions are related by multiplication on i. 

As in the two-magnon case we suggest that the reduced wave function should be 
bounded 

sup b (k, to, n) < oo. (69) 

m,n 

Despite the system (65), (67) is linear a proper generalization of the straightforward 
matrix approach used in the previous section is unclear for it. Instead one may treat (65) 
by the Fourier substitution 

bo(k, to, n) = <f(k, m)0(k, n), (70) 
which results in the following two-parametric exponential solution 

b (k,m,n) = e i (« im+ « n >, (71) 
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related to dispersion 

E(k : q 1: q 2 ) = 2[3J X + J 2 (cosk 1 + cos/c 2 + cos/c 3 )], (72) 



where 



fci = g-5i, ^2=g+gi-g2, A;3=g+g2- (73) 



Since 



e ifc/3(m+n+p) 6 (M - m,p-n) = S kim+k2n+k3P \ (74) 

one can naturally associate (71) with the wave function of a triple of magnons related to 
wave numbers k±, k 2 and k%. 

Instead of g"i i2 we shall mainly use the parameters 

k 2 -h _ q 2 h-k 2 _ q-y 

9i = — g — = gi ~ 2"' 92 = — 2 — = Q2 ~~ 2"' ( ^ 

considering them as generalizations of the parameter q from Eq. (53). The pair q 12 is 
more convenient for representation of pure scattering states related to real &i,2,3 with < 
k\ < k 2 < k 3 < 2n (a generalization of Eq. (2)) because the latter system of inequalities 
in terms of q± t2 has a very simple form. Namely < q± t2 < n and < q\ + q 2 < n. 
Nevertheless due to a rather compact representation (71) the parameters q\ j2 still will be 
remained in exponential factors. They also will be used for classification of states with 
complex wave numbers (Eq. (82)). 

When qi 2 are complex numbers one have to treat carefully the condition (69) and take 
into account that the energy (72) must be real. These conditions result in 

3(&)>0, J = 1,2, (76) 
( sin [k/3 — — sin [k/3 + $t.(qi — 52)] cosh Q(q 2 )^ sinh £s(qi) 

= (sin [k/3 + U(q 2 )[ - sin [k/3 + - q 2 )\ coshSfo)) sinh3(g 2 ). (77) 

The dispersion (72) is invariant under permutations of k\, k 2 and k% or equivalently 
under the following transformations of q = (q±, q 2 ) 

^i(q) = <?i - 52), ^2(q) = (-52, Qi - 52), w 3 (q) = (-q 2 , -ft) 

^(q) = (92 - ft, w 5 (q) = (q 2 - q~i, q 2 ). (78) 

In fact these formulas give a representation of the three-elements permutation group S3. 
It may be readily checked that all ujj are generated by lo\ and cj 5 . Namely, 

UJ 2 = U 5 ■ U>i, CJ 3 = U$ ■ liJ\ • U$ = LO\ ■ U$ ■ U>i, UJ4 = iO\ ■ (79) 
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The symmetry (78) allows to generalize the solution (71) and suggest the following 
ansatze (q = (q\, q 2 )), 



b Q (k, q, m, n) = Ax(k, q ) e ^ im+9 ~ 2n ) - A 2 (k, q ) e ^ im+ ^-^ n ) 

+ A 3 (k, q) e < (-«« m+ ^-®) n ) - A 4 (k, q ) e - i ^ m+ ^ 

+ A 5 (fc,q)e i(( « 2 - 9 " l)m - 9 " in) - A 6 (k,q)e l ^ 2 -^ m+ ^ n) . (80) 

For S(gi 2 ) = the expression (80) agree with (69) while Eq. (77) is satisfied iden- 
tically. However even when one of q\p has an imaginary part some of the amplitudes in 
(80) must turn to zero in order to ensure an agreement with (69). Besides according to 
(77) real and imaginary parts of q% t2 should be interdependent. More specifically let us 
divide the sector Q(q lf2 ) > on five subsectors 

Vi = [3f(g 1 )>0,3f(&) = 0], V 2 = [9f(?i) = 0,3f(g a )>0], 

V 3 = [0 < = V 4 = [0 3(g 2 )], 

V 5 = [0 < 3f(g 3 ) < (81) 

Let Ji will be corresponding to V% set of l-s for which there should be Ai(k, q) = 0. At the 
same time Qj will be the corresponding additional condition on g 12 following from (77). 
For each % we may gather a triple W; = [V»; Ji, G%\- A straightforward analysis based on 
Eqs. (77) and (80) results in the following classification 



Wi = 




3(g 2 ) = 0; {4,5,6}; K(g 2 ) = 2^)], 


VVi = 


> ; 


3f(g 2 ) = 0; {4,5,6}; &(q 2 ) = 2k/3 + tt], 


w 2 = 




3(g 2 )>0; {2,3,4}; = 2&(g 2 )], 


w 2 = 


[9f(?0 = 0, 


9f(g 2 ) > 0; {2, 3, 4}; = -2A;/3 + tt], 


w 3 = 


[0 < 


= 3(&); {3,4,5}; 3fc(&) = -&(&)], 


m = 


[0 < 9f(g0 


= {3, 4, 5}; = tofa) + tt - 2fc/3], 


w 4 = 


[0 < 3f(g x ) 


<3(&); {2,3,4,5}; Eq. (77)], 


m = 


[0 < 9f(g 2 ) 


<9f($i); {3,4,5,6}; Eq. (77)]. 



(82) 

Each Bethe state with complex g-s corresponds without fail to one of the W-s presented 
in (82). 

The system (65) (at m,n > 1), (67) exactly coincides with the well known one inherent 
in the XX Z model [13] . [14] . Nevertheless we shall give its solution within the ansatze 
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(80) in order to illustrate the straightforward approach used in the next section for S — 1 
and S = 2. 

Let us begin with pure scattering states for which the duality transformation (68) 
results in 

V(Mk, q )) = A l ^(k,ci), (83) 

where Ai(k,q) = Ai +e (k, q) for / = —2,-1,0. The 5 = 1 and S — 2 analogs of this 
formula are be used in Appendix for enumeration of three-magnon Bethe states. 
Substitution of (80) into (67) produces a linear system on the amplitudes Ai(k,q) 

6 

J2Mf ) (k, q )A l (k,c l ) = 0, (84) 
i=i 

where nonzero entries of the 6x6 matrix M^(k, q) are 



Mg>(*,q) = 


-Mif(k, q ) 


= Z^u^q),^), 


M 2 ( 2 0) (A;,q) = 


-M 5 ( 6 0) (A;,q) 


= Z(k, q.AO, 


M^(k,q) = 


-< } (M) 


= Zfau^AJ, 


M^(k,q) = 


-MH\k, q ) 


= Z(k,u 2 (q),A 1 ), 


M^(k, q ) = 


-M 2 ( 3 0) (A;,q) 


= Z(fc,W3(q),Ai), 


M$(k, q ) = 


-M 3 ( 4 0) (A;,q) 


= Z(fc,w 4 (q),Ai). 



(85) 



Here 



Z(ife, q, A) = cos + ^ gl ) - Ae^ 1+92 \ (86) 
while according to (75) and (78) 

^i(q) = (<?i + 92, -92), ^2(0) = (-91 - 92, 9i), w 3 (q) = (-92, -qi), 

^(q) = (92,-91-92), ^5(0) = (-9i,9i + 92)- (87) 

Since 

det M<°>(fc, q) = J] ^ - II M S+i^> ^ = °» ( 88 ) 

71=1. .6 n=1..6 

(here M$(k,q) = M^\k,q)) the matrix system (84) is solvable. Namely 

3 

i4 l (fc > q)=IjM I %_ i (fc > q), (89) 

7=1 

where M^\k, q) = M^+^A;, q) for Z = -2, -1, 0. 
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States with complex qi^ may be obtained from (89) by analytic continuation with re- 
gard to conditions presented in the list (82). It may be readily shown by a straightforward 
calculations that there are no solutions related to VVi,2,3 and W 4 ,5. This statement is a 
special (related to the three-magnon sector) confirmation of the string hypothesis proved 
for the XX Z chain [13], pi]. 

It may be readily proved that the S = 3 case is analogous to the S = one. It is only 
necessary to improve the representation (64) (in order to obtain the state with total spin 
S = 3) and replace everywhere A x on A 2 . 



5 S = 1 three-magnon sector 

General S = 1 three-magnon state has the following representation 



[3, 1, jfe) Xlm<n<p e 



,ifc(m+n+p)/3 



b^\k, n — m,p — 



n 



\l)^...\l)i...\l)l... + bf\k,n-m,p-n)...\l) b m ...\l) a n ...\l) 



+bi(k, n — m,p — 



n 



|i>! 



|i>^-.-|i>S- 



(90) 



and depends on the three wave functions b{~' (A;, to, n). At m,n > 1 the Schrodinger 
equation for bi' 3 \k,m,n) separates on three independent linear subsystems of the form 
(65) (one have only to replace bo{k,m,n) on bf' 2 ' 3 \k, to, n)). However for m, n— 1 one 
gets 

(6Ja + J 2 + p 3 ^\k, l,n) + J 4 bf\k, l,n) + J 2 ^e~ ik/3 b^\k, l,n-l) 
+e ik/s b a) i ; n + l) + e^/s^ 1 ) (k,2,n- 1) 
+e- ife / 3 &f ) (£;,2,n)l = Eb^\k, l,n), 
(6Ji + J 2 + |j 3 )&£ 2) (A;, l,n) + J 4 6? ) (fe, l,n) + J 2 [e~ ifc / 3 6? } (A:, l,n - 1) 

+e ik/3 bf ) (ife, 1, n + 1) + e ifc/3 ^ 2) (jfe, 2, n - 1) 
+e -ife/3 & (2)^ )2)n )l = Ebf\k,l,n), 

(4Ji + £ )&S 3) (A;, l,n) + (J 6 - ^X&f^, l,n) + &f(fc, l,n)) 
+ J 2 [e' ik/3 bf ] (jfe, 1, n - 1) + e^^bf ] (k, 1, n + 1) + e ife / 3 6S 3) (jfe, 2, n - 1) 
+e- ife / 3 fof ) (Jfe,2,n)] = £&S 3) (£;, l,n), 
(4 Jj + e )b? (k, to, 1) + ( J 5 - J 4 ) [&S 2) (k, to, 1) + 6< 3) (&, to, 1)] 
+J 2 [e ifc/3 6S 1) (A;,m - 1, 1) + e - **/ 3 ^ (ife, to + 1, 1) + e' ih/3 b ( i\k,m — 1,2) 
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+e ifc/s 6 (i) ( m> 2)] = Eb? (k,m, 1), 
(6Jx + J 2 + |j 3 )6S 2) (A;,m, 1) + J 4 bf\k,m, 1) + J 2 [e ifc / 3 &S 2) (A;, 



m 



1,1) 



+e- ik / 3 bf ] {k,m + 1,1) + e- ik ^bf ] {k,m - 1,2) 
+e ifc / 3 &S 2) (A;,m, 2)] = Ebf\k,m, 1), 




m 



+e -ifc/3 fe (3)^^ m+ i ? i) _|_ e - <fc / 3 6i 3) (A;, m — 1,2) 
+e ik / 3 bf\k,m, 2)] = Ebf\k, m, 1). 



(91) 



Introducing again the unphysical values ftp (k,m,0) and 6^ (/c, 0, n) one can reduce 
(91) to the form (65) by producing the following system of Bethe conditions 



(A 2 + A 1 )b?(k, l,n) + (A 2 - AJb^ik, l,n) = e ik / 3 b? (k, 0, n) 

+e~ ik/3 b ( ±\k,0,n + 1), 
(A 2 + AO&f^fc, l,n) + (A 2 - AO&^fc, l,n) = e lk / 3 bf\k, 0, n) 
+e- lfc / 3 6f } (A;,0,n + l), 
2Ac4 3) (A;, l,n) + §(A - A 2 )[b?(k, l,n) + bf\k, l,n)} 

= e ik / 3 bf\k,0,n) + e'^bf (k,0,n + 1), 
2A &[ 1) (A;,m, 1) + |(A - A 2 )[&S 2) (A;, m, 1) + bf\k,m, 1)] 
= e- ik / 3 bP(k,m,0) + e ik / 3 b^ (k,m + 1,0), 
(A 2 + A 1 )bf ) (k, m, 1) + (A 2 - A 1 )bf ) (k, m, 1) = e~ ik / 3 bf } (fc, m, 0) 

+e lfe / 3 6S 2) (A;,m + l,0), 
(A 2 + AJbf (k, m, 1) + (A 2 - Ai)6f ) (A;, m, 1) = e" ifc / 3 ^ 3) (A;, m, 0) 



As it was in the S = case the system (92) (as well as the three separate subsystems 
of the form (65) for bi \k,m,n)) is symmetric under a duality transformation 



Before developing a general analysis of the system (92) we shall at once find all the 
cases when it may be reduced to the XXZ-type form (67). 
First of all for 



+e ik / 3 bf\k,m + 1,0). 



(92) 



£>(&n*;,m,n)) 



b[ (k,n,m). 



(93) 



A, 



o — 



A 



i — 



(94) 
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the system (92) decouples into three XXZ-type subsystems (67) and therefore is com- 
pletely solvable. Namely 

bY'(k,q,m,n) = ajb (k, q,m,n), j = 1,2,3, (95) 

where 01,2,3 is a triple of arbitrary parameters. Labeling the coupling constants related 
to (94) by an upper index "(0)" one may readily obtain from (32) 

t(o) 

7(0 _ t(0 t(0) _ a t(0) a ._ Jrr , , 

J d - ~ J i 7 J u - 4J / > A o,i,2 - ~r~M)i w 
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or equivalently 



Jf = Jf = 0, a£> 2 = 1 + ^. 07) 

In this case according to (17) and (97) an interaction between excited triplet rungs is 
spin-independent. In more detail the relation between this model and XXZ chain was 
studied in [27] . 

Besides the complete separable case (94) there are also two configurations of A-s for 
which the system (92) possess a partial solution of the form (95) however with special val- 
ues of the ratios aii/aj = 1, 2, 3). Indeed substituting the ansatze (95) into the system 
(92) one readily makes sure that the latter may be reduced to (67) with an appropriate 
parameter A only under the following system of conditions 

(Ax + A 2 - 2A)«! + (A 2 - Ax)a 2 = 0, 

(A 2 - Ai)oti + (Ai + A 2 - 2A)a 2 = 0, (98) 

(Ax + A 2 - 2A)« 2 + (A 2 - Ai)a 3 = 0, 

(A 2 - A!)a 2 + (A! + A 2 -2A)a 3 = 0, (99) 

3(A - A)ai + (A„- A 2 )(« 2 + « 3 ) = 0, 

3(A - A)a 3 + (A - A 2 )( ai + a 2 ) = 0. (100) 

It may be easily observed that a trivial solution of the subsystem (98) may be nontriv- 
ially extended as a solution of the whole system (98)-(100) only in the case (94). From 
the other hand a nontrivial solution of (98), namely a.\ = a 2 , exists only for A = A 2 . 
Besides for A x = A 2 = A the system (98) is satisfied for all «i, 2 . Extension of these two 
solutions on the subsystems (99), (100) results in 

aii = «2 = «3, A = A = A 2 , (101) 
4ai = -a 2 = 4a 3 , A = Ai = A 2 . (102) 



(0) 
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Turning to the general (XXZ-irreducible) case we suggest the Bethe ansatze, 

$\k, m,n) = B[ j) (k, q) e i( -^ m+ ^ - B { j\k, q)<j^ m +^-&) n ) 
+ J B^ ) (^ ) q)e i( ~ 9 " 2m+{9 " 1 "^ 2)n) - J B^' ) (A;,q)e- ifem+ « in) 
+B^\k, q) e i (W2-gi)m-9in) _ B$\k, q)e i (^-^) m +^ n ) . (103) 

Classification of states with complex q 12 has the form (82). However each Ji in (82) 
is now a set of l-s for which all B^\k, q) = 0. In the present paper we shall not study 
S = 1 and S = 2 three-magnon Bethe states with complex wave numbers. For the pure 
scattering states the duality (93) reduces on the amplitudes as follows 



V(By\k, q )) = Bi j _\(k,a), 
Substitution of (103) into (92) gives 



(104) 
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J2M^(k,q L )B l (k,q)^0, 

i=i 

where the vector column Bi(k,q) for / = 1, 18 is defined as follows 

B 6U _ 1)+m (k,q) = B$(k,q), j = 1,2,3, m=l,...,6, 
while nonzero entries of the 18 x 18 matrix M^ l \k,q) are the following, 



(105) 



(106) 



M, 



M. 



(i) 
ii 

(i) 

22 



M. 



(i) 



33 



M, 



(i) 



44 



M, 



(i) 



55 



M, 



(i) 



(in 



M. 



(i) 



77 



M, 



(i) 



M- 



M- 



99 

(1) 
11,11 

(1) 
14,14 



fc,q 
k,q 
k,q 
k,q 
k,q 
k,q 
k,q 
k,q 

k,q 
k,q 



= -M$ 17 (k,q)=Z(k,u 5 (q),A Q ), 



Mg\k,q) = -M[l] 12 (k,q) = -Mg^(fc,q) = z(k,q, ^y^)' 



M[l] 13 (k,q) = Z(k,LU 1 (q),A ), 



Ai + A 2 



= M« 10 (A;,q) = -Mff^q) = -M« 14 (A;,q) = z(fc >W2 (q), 
= -M 1 « 15 (A;,q) = Z(A;,cu3(q),Ao), 
= M« 12 (A;,q) = -Mg (fc,q) = -M$ w (k,q) = z(k, u 4 (q), ^y^) 

)• 
)• 
)■ 



i< 13 (M) 
Mg 15 (A;,q) 



-M«(A;,q) 
-M«(A;,q) 



-M« n (A;,q) = z(A;,^ 5 (q), 
-M 1 «(A;,q) = zfA;^ 1 (q), 



Ai + A 2 
2 

Ai + A 2 



Mg; i7 (A;,q) = -M«(A;,q) = -Mg\k,q) = z(k,u 3 (q), 
-Mg\k,q) = Z(q,A ), 



Ai + A 2 



20 



M$(k,<C 

Mi%(k,q, 
M^(k,q) 

M^(k,q] 



M$(k : 



<io(^q; 



M (1) 
M (1) 



Mi%(k,q) 



q! 

q; 



= -M« 5 (A;,q; 
= -M« 4 (A;,q; 
= -Mg 1 (fc,q: 
= M 3 ( j 1 1 ) 6 (A;,q) = -M 1 ( j; 6 (A;,q; 
= M 5 ( j \ ) 7 (A;,q) = -M 1 ( i;3(A; ! q; 
= M 5 ( j 1 1 ) 8 (A;,q) = -M 1 { i; 2 (A;,q; 
= M«(A;,q) = -M 1 ( 1 i ; i8 (A;,q; 
= Mg\k, q ) 
= Mg^(fc,q) = 
= Mg; B (fc,q) = 
= M 1 «(A;,q) = 
Mg! 6 (fc,q) = 



q" 



= -M[\] 17 (k, 
-Mi%(k,q] 
-Mi%(k,q] 
-Mi%(k,q] 



>(k,q, 



A 2 - An • 
2 U P »<?2 



= -Mg } (A;,q) = Z(A;,cu 2 (q),Ao), 

= -M 3 ( p(A;,q) = Z(A;,^ 4 (q),Ao), 

= M&fo q) = -Mg^fc, q) = -M^ n (k, q) = ^ 

= M^ 4 (A;, q) = -M^ffc. n^l = -M, ( iL(fc. = A ° ~ A V ig2 

= <i } 5 (^q) 



■e ,w , 
■e" 



-< ) u (*,q) = ^ 



^ - A °~ A2 ^ 



-Mg; i2 (fc,q) = ^5_ 

A 



■iqi 



= -M«(A; ! q) = ^ y ^e-^^), 
= -M«(Jb,q) = ^^e^), 

= -M$ u (k,q) = A l_A e -^i + , 2 ) ) 
= -Mg; g (fc,q) 



\ A 2 -Ai 



Ai - A 2 
2 

A 2 - Ai ,„ 



, m „ * Ai-A 2 



■e i<Zl , 



•Mg 10 (A;,q) 2 x , 

of the 5 = 1 problem it is necessary to 
. (92). In the Bethe Ansatze framework 

(108) 



2 



As in the case (94) for complete solvability of the S = 1 proble 
have three independent solutions of the system (92). In the Bethe 
(103), (105) this results in 

rank(M (1) (A;,q)) = 15, 



-e~ iqi . 



(107) 



and therefore in 

P«(A;,q)=0, n = 0,l,2. 

Here Pi 1 \k,q) are coefficients of the characteristic polynom 

18 

c,q)-A/| = ^Pi 1 )(A;,q)A". 

n=0 



Direct 



\M w (k, 

n=u 

MAPLE gives 



calculation based on the computer algebra system 
P (1) (A;, q ) = detM (1) (A;, q) = 0, 



(109) 



(110) 



(111) 
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therefore even in the general case rank(M^ (k, q)) < 17 and the system (105) always has 
at the minimum one solution. Its general form is represented in the Appendix A. 

For the next coefficient P^ 1 \k,q) we have obtained by machinery calculations the 
following factorization 

P[ l) (k, q) = A (i _ e 2 ^) 2 (l - e 2 ^) 2 (l - e 2 ^+^) 2 p?\k, q), (112) 

where 

P«(A;,q) = e- l ( llfc+3l91+31(?2 )/ 3 ^ Q m ,„, p (A , Aj, A 2 )e t(mk+n(ll+pc >^ 3 . (113) 

m,n,p>0 

The sum in (113) contains 95052 terms (that is why P^\k, q) can not be represented in 
the format of this paper). According to (112) and (75) the condition (109) is satisfied at 
n — 1 if either any two wave numbers from hi, k 2 and k 3 coincides or at 

P 1 (1) (A;, q )=0. (114) 

The former three cases are similar to the case k\ = k 2 in the two-magnon problem studied 
in the Section 3. Three-magnon solutions of this type will not studied in the present paper. 
Turning to the Eq. (114) we shall confine ourselves by the problem of its solvability for 
all wave numbers. Namely we shall postulate 

Qm,n,p 

(A ,A 1 ,A 2 ) = 0, (115) 

to be valid at all m,n,p. 

Despite the system (115) depends only on a triple of unknown variables it is practically 
unsolvable by the MAPLE Grobner package on a personal computer with RAM about 2 
Gb. Luckily (as it may be directly checked by machinery calculation) 

2g8 ' 4 ' 13 25 92A 2 A 2 Ai g4 ' 8 ' 15 = 12(A ° " Al) ' + 15(Al " ^ + ~ Ao) ^ (U6) 

so except (94) there are no solutions with A AiA 2 ^ 0. 

In each of the three cases A ,i,2 = the reduced system (115) is essentially simpler 
than the initial one and may be readily solved on the personal computer. Calculations 
based on the Grobner package gave four pairs of solutions. We shall represent them as 
sets of A-parameters: A = [A ,Ai,A 2 ], and additionally as the corresponding sets of 
the coupling constants: J = [Ji, Jd, J rr , Jn, Jdd] and J = [J 2 , J3, J4, J5]. Note that the 
parameters J r and J\ remain indefinite. This is rather evident because both of them 
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correspond to the term proportional to Q in the Hamiltonian. But according to (23) the 
former has no affect on the Bethe equations. Namely the solutions are the following 



A^) = [±1,0, ±1], J^ = [1,0,0,4,0], = [±1,0,1,1], (117) 

A^) = [0,±1,0], j( 2 ' + ) = [1,0,-2,4,0], J^) = [±1,-2,1,1], (118) 
A( 3 '±) = [0, ±|, ±|], J< 3 > + ) = [1,0,4,0,-4], JP.±) = [±1,4,0,-1], (119) 

A (4 < ±} = [±~,0,0], J< 4 '+) = [1,0,1,0,-4], ^ = [±1,1,0,-1]. (120) 

It may be readily shown that the condition (27) is satisfied only for " +" solutions, while the 
"-" ones may be obtained from them by the symmetry (28). That is why we have omitted 
representations for j( 1 ' 2 ' 3 ' 4 ' - ). However they may be readily obtained from j( 1 ' 2 ' 3 ' 4 ' - ) using 
Eqs. (20). 

The models related to A 1,+ and A 2 ' + were first presented in [28]. Then the former 
one was intensively studied in [3]. Algebraic structures related to A x ' + and A 3 '~ models 
as well as to the model (94) with A 0j i )2 = 1 were presented in [23]. (However the cases 
(118), (120) and the general case (94) were not discussed in [23]). 

As it follows from (101) and (102) all the models (117)-(120) have a XXZ-type solution 
(95). The remaining pair of solutions may be chosen in different ways. (In other words we 
do not know the best choice of basis in the two-dimensional solution subspace additional 
to (95)). The basises obtained by machinery calculations within MAPLE are presented 
in the Appendix. 



6 S = 2 three-magnon sector 

A S = 2 three magnon state related to S z = 2 has the following form 



|3, 2, A:)" = J2 

m<n<p 



e ik(m+n+p)/3 r 



b^(k,n - m,p - n) . . . |1)+ . . . (|1>+ . . . |1) 3 - |1) 3 . . . |l>+)„. 
+ ^(fc,«-Ti»,p-n)...(|l>+...Il>» -|1>^...|1>+)...|1>+... (121) 

For m, n > 1 the amplitudes b 2 ' (k, m, n) separately satisfy the Schrodinger equation 
(65) while for m = 1 or n = 1 

(4Ji + e 2 )b ( 2 ) {k, l,n) + J 2 [e- ik / 3 b£\k, l,n- 1) + e ik ^b { 2 l \k, 1, n + 1) 
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+e ik/3 b m ^ 2 , n - 1) + e- ifc / 3 4 1} (k, 2, n)} = Eb 2 1] (Jfe, 1, n), 
(4Ji + £i)4 2) (^, l,n) + J 4 fc? ) (A:, l,n) + J 2 [e" ifc / 3 4 2) (A;, l,n - 1) 
+e ifc/3 ^ 2) (Jfe, 1, n + 1) + e ik/3 b { 2 2) (k,2,n - 1) 

+ e - ifc / 3 ^ 2) (M,n)] = M 2) (^> 1,"), 
(4Ji + e 2 )b 2 2 \k,m,l) + J 2 [e ik / 3 b {2) {k, m - 1, 1) + e- ik / 3 b { 2 2) (k, m + 1, 1) 

+e _ ifc/ 3 6 (2) ^ m _ ^ 2 ) + e *fc/3 6 (2) m? 2 )] = ^(2) ^ m? ^ 

(4Ji + £i)^ 1} (A;,m, 1) + J 4 b { 2 \k,m, 1) + J 2 [e ik/3 b 2 \k,m ~ 1, 1) 
+e -ifc/3 6 (i) m + 1? !) + e -ik/z b W (k,m- 1,2) 

+e ik/3 b 2 1 \k,m, 2)] = Eb 2 1 \k,m, 1). (122) 

Introducing again the unphysical amplitudes we obtain from (122) the corresponding 
system of coupled Bethe conditions 

2A 1 4 1) (Jfe, m, 1) + (A 2 - Ai)4 2) (A;, m, 1) = e'^b? (k, m, 0) + e ik / 3 b { 2 1] (k,m + 1,0), 
2 A 2 4 2) (Jfe, m, 1) = e" ifc / 3 6f (/c,m, 0) + e ik / 3 b {2) (k,m + 1,0), 
2Ai4 2) (Jfe, 1, n) + (A 2 - A^ (Jfe, 1, n) = e lk / 3 b {2) (Jfe, 0, n) + e~ ik / 3 b 2 2) (Jfe, 0, n + 1), 

2A 2 4 1} (A;,l,n) = e ik / 3 b 2 1] (k, 0, n) + e- ik / 3 b 2 1] (k, 0, n + 1), (123) 

invariant under duality transformation 

V(bf(k,m,n)) = b ( 2 ~ j \k,n,m). (124) 

For 

Ai = A 2 , (125) 

or (according to (32)) 

Ju = -4J d , (126) 

this system decouples into a pair of the XXZ-type subsystems (67) on b 2 1 ' 2 \k,m,n). In 
this case the general solution 

b 2 \k, q, m, n) — (3jb (k, q, m, n), j = 1,2, (127) 

depends on q and two arbitrary parameters Pip- 
it may be readily proved that the XXZ-type solutions (127) exist only under the 
condition (125). 
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In the general case making the standard substitution 

b { i\k,m,n) = C[ j \k, q )e i ^ m +^ n '> -C ( 2 j \k, qy^^+tfi-w)' 
+C^\k,q)e i ^ 2m+ ^ 1 -^ - C^\k, q) e -^ 2m +« in ) 
+C ( 5 j \k, q)e i (W2-gi)m-gm) _ C$\k, c^e i{ ^ 2 ~^ m+ ^ n \ 

one results in a linear system 

12 

J2M^(k, q )C J (k,c l )=0, 

3=1 

where as in (106) 

C 6 y_i)+ m (fc,q) = C<J } (A;,q), j = 1,2, m = l,...,6, 
and the 12 x 12 matrix M^ 2 \k, q) has the following nonzero entries 



<> 


[k, q) 


= 


-Mg; u (fc,q) = 


= Z(fc,w 5 (q),Ai), 


M 2 ( ? 


(fc, q) 




fOl 

-Mg; i2 (fc,q) = 


= Z(fc, q,A 2 ), 


M$ 


(fc, q) 




-M 1 ( ? 7 (A; ) q) = 


^(^(q^Ax), 


Mi? 


(fc, q) 




-Mg ) (fc,q) = 


Z(fc,w 2 (q), A 2 ), 


M 5 ( ? 


(A;, q) 




-M 8 ( ?(A;,q) = 


Z(fc,w 3 (q),Ai), 


M 6 ( ? 


[k, q) 




-M^(fc,q) = 


: Z(fc,w 4 (q), A 2 ), 


M<? 


& q) 




-M<?(fc,q) = 


Z(fc,w 5 (q),A 2 ), 


Mi? 


(fc, q) 




-M 5 ( ?(A;,q) = 




Mi? 


[k, q) 




-M 6 ( ?(A;,q) = 


Z(fc,Wi(q),A 2 ), 


M (2) 

JW 10,10 


(fc, q) 




-M 1 ( ?(A;,q) = 


Z(fc,w 2 (q), Ai), 


M (2) 
jw ii,ii 


% q) 




-M 2 ( ?(A;,q) = 


Z(fc,w 3 (q), A 2 ), 


M (2) 

Jw 12,12 


& q) 




-M 3 ( ?(A;,q) = 


Z(fc,w 4 (q),Ai), 


M{? 


(fc, q) 




-M 1 ( ? 5 (A; ! q) = 


Ai - Ao • 
2 


MS 


% q) 




-M 1 ( ? 4 (A;,q) = 


Ao - Ai 
2 


M (2) 

JW 39 






-M 1 ( ? 1 (A;,q) = 


Ai - Ao • 
2 


M (2) 

JW 3,10 


% q) 




-Mg; 6 (fc,q) = 


. A 2-A x 
2 
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r (2) 
5,11 

r(2) 



r(2) ( 

83 1 
r(2) , 



K^q) = -^(*,q) = ^e^) 



'1311 



As in the S = 1 case we shall concern only on the pure scattering states for which the 
duality (124) gives 

V(C®{k,q)) = C®{k,q). (132) 

To be completely solvable the system (129) must posses two independent solutions. 
Equivalently there should be 



rank(M (2) (fc,q)) = 10. 



(133) 



According to machinery calculation 

detM (2) (fc, q) = -6(Ai - A 2 ) 2 (l - e 2 ^ 1+92) ) 2 (l - e - 2iqi ) 2 (l - e - 2iq2 ) 2 Y 2 (k, q), (134) 



where 



Y(k,q) = (A x -A 2 ) 2 -1 



cos k + 



(A 1 + A 2 ) 2 -l 



cos ■ 



k - Aq± - 2g 2 



k + 2gi + 4g 2 fc + 2g x - 2g 2 
+ cos h cos 



-4AiA 2 (Ai + A 2 ). 



(135) 



A condition 



detM (2) (A;,q) = 0, (136) 
will be satisfied at all k, q± and g 2 either in the case (125) or in the four additional ones 

A x = ±l, A 2 = (137) 



and 



A 1 = 0, A 2 = ±l. 



(138) 



Machinery calculations show that in all these cases the condition (133) is satisfied. The 
corresponding solutions of the system (129) are presented in the Appendix. 



7 Integrability and the Reshetikhin condition 

A well known alternative to the Coordinate Bethe Ansatze is the so called Algebraic Bethe 
Ansatze or the Inverse Scattering Method [15]- p2]. It is based on the representation of 
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the finite dimensional matrix H related to the local Hamiltonian density H n>n+ i as a 
derivative of the corresponding i?-matrix. 

H = Atf(A)| A=0 . (139) 

The latter satisfies the Yang-Baxter equation, 

Rn(X - /i) J R 23 (A) J R 12 (/i) = RMR 12 (X)R 23 (X - /i), (140) 

and the initial condition, 

R(0) oc /, (141) 

(where again I is an identity matrix). 

From (139)-(141) follows the Reshetikhin condition [T7] . 

[H 12 + H 23 , [H l2 , H 23 \] = K 23 - K 12 , (142) 

which for the Hamiltonian density (17) with Jq = gives the following system of equations, 

J 2 J 4 (Ji + Js + J 5 ) = 0, 

J2(J4-J 5 )(J4 + h) = 0, 

J 2 (J 4 - J 5 )(2Ji + 2J 3 - J 4 + 5J 5 ) = 0, 

(J 4 - J 5 )(J|- J 5 2 + 2J 4 J 5 ) = 0, 

J 5 (J 2 2 -2J 4 2 - J 5 2 + 2J 4 J 5 ) = 0, 

J 2 (J 2 + 2J 1 J 3 -4J 2 + 4J 4 J 5 ) = 0. (143) 

Taking at the first J 2 = one gets from (143) J5 = 0. This case with degenerate one- 
magnon dispersion is of poor physical interest and was already studied in [26J. Taking 
now J 2 7^ and using (32) one can subdivide the system (143) on two subsystems, 

(A 2 - A )(A 2 -A 1 )(3A 1 -2A 2 - A ) = 0, 
(A 2 - A )[9(Aa- A 2 ) 2 -4(A 2 - A ) 2 + 9] = 0, 
(3A 1 -A 2 -2A )[9(A 1 -A 2 ) 2 + 4(A 2 -A ) 2 -9] = 0, (144) 



and 



(Ai - A 2 )(3A r + 3Ai + A 2 - A ) = 0, 
(A 2 -A )(3A r + A 2 -3A 1 + 5A ) = 0, 
3A r (3Ax + A 2 - A ) + 9A 2 -18A 1 A 2 + A 2 + 4A 2 A -5A 2 = 0, (145) 
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where A r = J r /(J; — Jd)- 

We have separated (145) from (144) because the variable J r related to the term propor- 
tional to Q is auxiliary and according to (23) has no affect on integrability. Nevertheless 
solvability of the Yang-Baxter equation (140) for a given H depends on the value of J r . 
Therefore however it is necessary for the system (144), (145) to be solvable a concrete 
value of J r obtained from it has no affect on integrability. 

The subsystem (144) has three solutions. The first one is the solution (94) for which the 
subsystem (145) is also solvable. The remaining two solutions of (144) are the following, 



A?) 2 



A = A 2 , (Aj 



Ai = A 2 , 4(A -A 2 ) 



9. 



(146) 
(147) 



A substitution of (146) into (145) shows that the latter subsystem is solvable with respect 
to A r only for, 

AiA 2 = 0. (148) 

Together (146) and (148) result in (117) and (118). 
Analogously a substitution of (147) into (145) gives, 



AxAo = 0. 



(149) 



Together (147) and (149) results in (119) and (120). 

Corresponding to the integrable cases i?-matrices were already presented in [21] within 
the following basis in the space C 16 



+3 



£j, fg+i = |0) ® ej, fu+i 



|0>, /i6 = |0>®|0). 



(150) 



Here i,j = 1,2,3 and t\ = |1} +1 , e 2 = |1)° and e 3 = |1) l . 

In this basis the i?-matrix corresponding to (94) has the block XAZ-type form 



#°)(A) 



sinh(A + r/)/ 9 

sinh rjl^ sinh A/3 

sinh A/3 sinh 77/3 



\ 



V 











(151) 



sinh (A + 77) 1 



For a very special value of r] it was also presented in 

In the cases (117) (for J r = 0) and (118) (J r = J{) the matrices H are correspondingly 
the normal and graded C 4 ® C 4 permutators V4 and V\. (In the latter case the subspace 
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generated by |0) has the negative grading). The related /^-matrices have a rather simple 
form, 

.R (1 ' 2) (A) = i]Ii e + XH. (152) 

Integrability of these models was first noted in [28J. The case (117) was intensively studied 
in®. 

The i?-matrices related to (119) (for J r = Ji) and (120) (2 J r = 5/;) also have block 
forms, 



( 



\ 



r{X,rj ) 












sinh r] I 3 sinh A/3 

sinh A/3 sinh 770/3 





sinh 7/9/3 sinh A/3 

sinh A/3 sinh r] I 3 











where sinhr/o = v5/2 and, 



r{\,r)o) 



N 




sinh (A + r/o) 1 





sinh (r/o - A) J 
\ 



f 





























/ 





























f-9 





9 





-9 

















/ 























9 





f-9 





9 























/ 

















-9 





9 





f-9 





























/ 





























/ 



(153) 



(154) 



V 

(/ = sinh (A + r/o), g = sinh A) . 

The matrix r(A,r/o) itself satisfies the Yang-Baxter equation and describes the S = 1 
biquadratic spin chain. As it was shown in [29] this i?-matrix as well as its generalization 
(related to arbitrary 77) are related to the Temperley-Lieb algebra. 
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8 Action of the S3 group on the eigenspaces 

As it will be shown below (see Eq. (162)) the S^-action (87) in the g-space results in 
corresponding symmetry of Bethe wave functions. The latter is useful (see Appendix A) 
for compact representation of amplitudes. 

First of all let us consider the case S = (which is analogous to S — 3). The matrix 
M^(k,q) possess the following symmetry 

M(°)(A^(q)) = 4%,0M(°>foq)jgV;), (155) 
where the matrices and J^\ujj) give left and right representations of the group 

Jf Jf = Jf (u* • 4\^)4\^) = J$\ Uj ■ Ui ). (156) 

Explicit expressions for the matrices J^R^j) ma Y be obtained from Eqs. (79), (156) and 
the following representations for generators 

\ o 5l i h J \ o 1>5 i / 

= - [ /2 ° 2 ' 4 I , 4 0) K) = ~h- (157) 

\ 4 ,2 /4 / 

Here by O mi „ we denote a m x n matrix with all zero entries while by I n a n x n matrix 
with units in the second diagonal (and all other entries equal to zero). 
Similar relations 

M^\k, Uj (q)) = J?' 2 \u j )M^\k,c l )J R 1 ' 2 \u; j ), (158) 

with 

4 1 i = ^s®4°j i , J^ R = h®J { ^ w (159) 

are also valid for M^ 2 \k,q) given by (107) and (131). 

The symmetry (158) allows to produce new solutions of the Eqs. (105) or (129) from 
the known one (for Eq. (84) the result is trivial). Indeed if 

M (n) (A;,q)F (n) (A;,q) = 0, (160) 

for some vector F^ n \k,q) (dim(F^(k, q)) = 18, dim(F( 2 )(A;, q)) = 12) then according to 
(158) 

MW(^q)jW(a; 3 )FW(fc,^(q)) = 0. (161) 
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In other words we have obtained the following action of the group S3 on the eigenspaces 

^ (fl) )(*,q) = 4 n) (^ (n) (^-(q))- (162) 

Here Uj{F^) is vector related to the new solution (which in fact may coincide with the 
present one). 

9 Summary 

In the present paper we analyzed two- and three-magnon problems for a rung-dimerized 
spin ladder. It was shown that the Bethe form of the two-magnon solution may be 
obtained in a straightforward manner from the corresponding Shrodinger equation. 

The three-magnon problem was first analyzed in general outlook in all sectors of total 
spin 5 = 0, 1,2,3. It was shown that at all S the reduced to the center of mass frame 
Shrodinger equation is invariant under the corresponding duality transformation while 
the Fourier substitution (70) naturally results in Bethe form of wave function. 

Applicability of the Bethe Ansatze for the three-magnon problem was analyzed sepa- 
rately in all sectors of total spin. It was shown that for S = and S = 3 the problem is 
always solvable and the corresponding solution has form typical to the XX Z model. The 
sector S — 1 is completely solvable in the five cases (94) and (117)-(120). Nevertheless 
a special partial solution (see Appendix A) exist for all values of the coupling constants. 
The sector S = 2 is solvable under one of the conditions (125), (137) or (138). Explicit 
expressions for the solutions are presented in the Appendix. 

The result was compared with the previous consideration based on an analysis of 
solvability of the Yang-Baxter equation. It was shown that the three-magnon problem for 
a Hamiltonian H is completely solvable within the Coordinate Bethe Ansatze if and only 
if the corresponding R- matrix exist for some Hamiltonian in the orbit H + aQ (a is real). 

Finitely it is shown that the S'3-symmetry of the Bethe Ansatze equations results in 
the action (162) of the group S3 on the space of Bethe vectors. 

The author is very grateful to P. P. Kulish and M. I. Vyasovsky for helpful discussions. 

A Partial solution in the S = 1 sector 

An explicit form of the special partial solution of Eq. (105) obtained by MAPLE is rather 
complicated. For example the expressions for Bj(k,q) at j = 1,..,6 and j = 13, ..,18 
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contain 1106 terms while the expression for Bj(k, q) at j = 7, .., 12 contain 1090. 

Since this solution is in general a single one it must be S^-symmetric and auto- (or 
anti-auto) dual. It may be readily proved that these symmetry properties allow to obtain 
all components from Bi(k,q) and B 7 (k,q) using Eqs. (104) and (162). Below we give 
representations for these two components. 

First of all Bi(k, q) possess the following decomposition 

B^q) = B[ s \k,q) + B[ a \k,q), (163) 
where the term B[ s \k, q) is symmetric under the transposition, 

k -> -k, qi <-> q 2 , (164) 

while B[ a \k, q) is antisymmetric. 

For B[ s \k, q) we have the following representation 

B[ s \k,q) = yF(A 1 ,A 2 ,A ) + yF(A 2 ,Ao,A 1 ) + ^(A ,A 1 ,A 2 ) 



Wiuiu 2 uz + W 2 h W 3 1 — — H )u 3 

-^1 ^2 ^1 ^2 ^1 ^2 " 



+ UiU 2 Uz 

+ ^ + ^(M)+S]. 



where 



F(A,A',A") = (A / -A)Z(fc,q,A w )Z(fc jW2 (q),A // )Z(fc,a;4(q),A w ) 

• [wim 2 «3 + AA' cos fc + AA'(A + A')]. (166) 

The parameters 

fk qi\ (k q 2 \ 

m = cos (- + -), u 2 = cos 

us = cos(^ + ^), z 3 =^l\ (167) 

have more simple form being expressed from q lj2 . 
The coefficients Wj for j = 1,2,3 are the following 

W 1 = 27A? - 5A^ + 8Ajj + 45A?A 2 - 75AxA2 - 90A2A 

+ 60A 2 A^ - 18A2A - 12AiAg + eOAxAsAo, 

W 2 = 18Af A 2 - 10Ai A^ - 45A?A - 50AiAg + 15A^A + 42A 2 Ajj 

+ 30(3A^ - Af)Al + 3A!A2A (65A2 - 39A! - 36A ), 

W 3 = ^2 + 15(A 1 -A 2 )(A 1 -Ao)(A2-A )(4Ao-3A 1 - A 2 ). (168) 
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For j = 4, 5, 6 they may be obtained from (168) by the following formulas (observed purely 
empirically) 

<p{W 2 ) w <p(W 3 ) 



Wa A!A 2 A ' W " Ai^Ao' 
where the homomorphism tp is defined as follows, 

AiA 2 A 



W 6 = (piWj, 



For B[ a \k,q) we found the following representation 



B ( i\K q) = f (A: - A 2 )(A X - A )(A 2 - A )u 3 B[ a >(k, q) 



15 



>(<*), 



where 



£j a) (£;,q) = MlM2 [(3-4A 1 A 2 -6A A 2 -2A A 1 )(^-^|) 



+ 3isinA; + 3i(l + 4(AoA 1 + A 1 A 2 + AoA 2 )) 



V3_ 

ZlZ 2 



2*(2A + 3A 1 + A 2 )(^ 
V z 1 



UlZ 2 , M2^l\ 



+ 



^T 3 



2i(A + 2A 2 ) 



■UiW 2 + u 2 v 1 



Z\Z 2 



UiV 2 + U 2 V\ 



2^2 



Z7Z. 



1*2 



and 



Vl 



Sm| 3 + 2 



^2 



'fc g 2 

sin 

\3 2 



6iA AiA 2 

'A; q 2 ~Qi 



V3 



sm (- + 



Representation of B 7 (k, q) is similar to (165) 

B 7 (k,q) = 45F(A ,A 2 ,A 1 ) + 27F(A ,A 1 ,A 2 ) 



+ U\U 2 U 3 



WiUiu 2 u 3 + W 2 - 



~ UiU 2 ~ /UiZ 2 M 2 ZA 
Vo 1- H/q I = 1 — I 



ZlZ 2 



+ 



4-r 



~ U 3 ~ (Ui u 2 
Z\ Z 2 V Z-i z 2 



+ 



We 

z\z\ 



where 



W 1 



Wo 



W 3 



(169) 



(170) 



(171) 



(172) 
(173) 



(174) 



= 27A? + 5A^ -32A^ + 45AiA 2 (A 2 - A x ) 

+ 72A 1 A (A 1 -Ao) + 120A 2 A (A 2 -A ), 

= 45A?A - 72A^A 2 - SOA^Ai + 75A 2 ! A + 20A^Ai + 12A^A 2 

+ 60(2A^ - A^Ag - A^A 2 ) + 3A 1 A 2 A (104A - 29Ax - 75A 2 ), 

= W 2 + 90(A 1 -A 2 ) 2 (A 1 -A )(A 2 -A ). (175) 



Again the parameters W^^ may be obtained from W\^,z according to (169) and (170). 
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B Additional 5 = 1 solutions 



We shall use here the following notations 



mi(fc,q,A) = Z(k,u 5 (q L ), A), m 2 (k, q, A) = Z(k, q, A), 
m 3 (k, q, A) = = Z(k, wi(q), A), m 4 (A;, q, A) = w 2 (q), A), 
m 5 (k,q,A) = Z(k,u 3 (q), A), m 6 (fc, q, A) = Z(fc, w 4 (q), A), 



(176) 



(for definition of Z(k, q, A) and cjj(q) see (86) and (87)). 

For A = A 2 = 1, Ai = the space of additional to (95), (101) solutions is generated 
by the vector 



B[ 1] (k,q 

B { 2 \k,q 
q 

B±\k, q 
q 

5«(A;,q 

i?; 2) (A;,q 
B^(k,q 

B?\k,q 
^ 3) (fc,q 
5, 0) (fc,q 



A; + 2g 1 + 4g 2 . fc - 4 ?1 - 2q 2 

2m A (k, q, l)m 6 (fc, q, 1) sin sin . 

6 6 

n u n . k + 2 qi + Aq 2 . k-A qi -2q 2 
2m 1 {k, q, l)m e (k, q, 1) sm sin . 

-m^k, q, l)m 2 (k, q, l)m 6 (k, q, 1), 

— mi(fc, q, l)m 2 (/c, q, l)m 3 (/c, q, 1), 

2m 2 (k, q, l)m 3 (A;, q, 1) sin ^ — sin ■ 1 



6 



6 



n u n u . k + 2q x + Aq 2 . k + 2 qi -2q 2 
2m 3 (k, q, l)m 5 (k, q, 1) sm : sm ■ 



6 

. . , , v • /c + 2g!+4g 2 . k + 2 qi 
2% sm (gi + q 2 ) sm sm 



2q 2 



6 6 
-im^k, q, l)m 6 (A;, q, 1) sin ( qi + q 2 ), 

-im 2 (k,q, l)m 3 (k,q, 1) sing 2 , 

. n . . fc + 2gi + 4g 2 . fc - 4 qi - 2q 2 
2im 3 (k, q, 1) sm q 2 sm sm . 

m 6 (k,q, l)sing 2 sin (q x + q 2 ), 
m 3 (k, q, 1) sing 2 sin (q 1 + q 2 ), 
0, (j,/) = (2,3-4), (3,2-5), 



m 6 (fc,q,l), 



(177) 



and its dual. 

For A = A 2 = 0, Ai = 1 the space of additional to (95), (101) solutions is generated 
by the vector 



^ 1} (M) 



B^\k, q) = m 5 (k,q, l)m 6 (fc,q, 1), 
(k, q) = -2m 6 (k, q, 1) sin + 2 ^ 



4<? 2 . k — Aqi - 2q 2 
sm . 
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(k,q) 


Bf 


(k,q) 


Bf 


(*,q) 


r(2) 


(«, q) 


B f 




Bf' 


(*,q) 




(*,q) 



/c + 2g!-2g 2 . fc - 4g x - 2g 2 
-2m 5 {k, q, 1) sm sm ■ 



6 



= im 6 (k,q,l)sm(q 1 + q 2 ), 



. . , % . fc + 2<?i - 2g 2 . fc - 4gi - 2g 2 
-2? sm (gi + g 2 J sm sm 

. k + 2gi + 4g 2 . fc - 4g x - 2g 2 

—22 sm g x sm sm , 

6 6 

im 5 (k,q, 1) sin gi, 

Sj 3) (/c, q) = -singi sin (gi + g 2 ), 

0, CM) = (2, 1-2), (3,1-3), (3,6), 



(178) 



and its dual. 

For Ai = A 2 = 3/2, A = the space of additional to (95), (102) solutions is generated 
by the vector 



B. 
B- 



(i 



(i 



B \ 



(2 



B. 

B 
B 
B 
B 



si 



(3 



B. 



(3 



B. 



(3 



5 



(3 



B, 



k,q 
fc,q 
/c, q 
fc,q 
fc,q 
fc,q 
fc,q 
fc,q 

fc,q 

fc,q 

fc,q 

fc,q 
k,q 



= -2m 5 (k, q, 3/2) singi sing 2 , 

= —2m 2 (k, q, 3/2) singi sing 2 , 

= m 4 (fc, q, 3/2)m 5 (fc, q, 3/2)m 6 (fc, q, 3/2), 

= mi(/c, q, 3/2)m 5 (fc, q, 3/2)m 6 (fc, q, 3/2), 

= mi(fc, q, 3/2)m 2 (£;, q, 3/2)m 6 (£;, q, 3/2), 

= mi(fc, q, 3/2)m 2 (fc, q, 3/2)m 3 (fc, q, 3/2), 

= -Bf(k,q) = m 2 (k,q,3/2)m 3 (k,q,3/2)m 4 (k,q,3/2), 

= -Bf (k,q) = m 3 (k,q,3/2)m 4 (k,q,3/2)m 5 (k,q,3/2), 

'k-q-i.- 2q 2 



-m A (k, q, 3/2)m 5 (fc, q, 3/2) 
-mi (A;, q, 3/2)m 5 (k, q, 3/2) 
-m 2 (fc, q, 3/2)m 6 (£;, q, 3/2) 



cos 



cos 



cos 



k-qi- 2g 2 



fc + 2gi + g 2 \ 
3 / 



— cos gi 

— cos gi 

— cos g 2 



^qi 



-IQ2 



-m 2 (k, q, 3/2)m 3 (£;, q, 3/2) cos y ^ J - cos g 2 

0, (j,/) = (1,1), (1,4-6), 



(179) 



-3/2, Ai = A 2 = the space of additional to (95), (102) solutions is 



and its dual. 

For A = 
generated by the vector 

B[ 1 \k,q) = Bf(k,q) — -im 5 (k,q,-3/2)m 6 (k,q,-3/2)smq 2 , 
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B^(k,q) = Bi 1 \k,q) = im 1 (k,< l ,-3/2)smq 1 

( k + q 2 - qi 
• [cos + cos(ft + ft) + - J, 

BP(k,q) = mi(A;,q,-3/2)m 5 (A;,q,-3/2)m 6 (A;,q,-3/2), Z = 1,2,3,4,5,6, 
B { *\k, q) = Bf\k,ci) = imi(fc,q, -3/2)m 6 (fc, q, -3/2) sin (ft + ft), 
Sf(fc,q) = (fc,q) = im 5 (A;,q,-3/2)sinft 



Bi j \k, q) = 0, (j,/) = (1,2-3), (3,1-2), (180) 
and its dual. 



A; + 2ft + g 2 e^ 2 \ 
cos + cosft + 



C S = 2 solutions 

For A x = ±1, A 2 = the space of solutions is spanned on 



dp 


(*,q) 


= q) 


= -im 5 (k,q)m 6 (k, q, ±1) sin ft, 




(*,q) 


= q) 


= 0, 


dp 


(*,q) 


= ^(^q) 


= im,i(k, q, ±l)m 2 (/c, q, ±1) shift, 


dp 


(*,q) 


= Cffoq) 


= ±m 1 (k, q, ±l)m 5 (fc, q, ±l)m 6 (fc, q, ±1), 


cf 


(*,q) 


= Cf(*,q) 


= ±mi (A;, q, ±l)m 2 (A;, q, ±l)m 6 (/c, q, ±1), 




(*,q) 


= ^ 2) (A;, q ) 


= ±m 2 (k, q, ±l)m 5 (k, q, ±1) (V^ 291 " 2 ^/ 6 



=p e^ 291 - 2 ^ 6 ), (181) 

and its dual. 

For Ai = 0, A 2 = ±1 the space of solutions is spanned on 



dp 




= dP(k,q) = 0, 


dp 




= im 5 (k, q, ±l)m 6 (A;, q, ±1) sinft, 


dp 


%q) 


= im 2 (k : q, ±l)m 6 (A;, q, ±1) sinft, 


dp 


[k,q) 


= im^k, q, ±l)m 6 (k, q, ±1) sin (ft + ft), 


<$> 




= im 4 (k, q, ±l)m 6 (fc, q, ±1) sin (ft + ft), 


dp 




= dP{k, q, ±1) = =Fm 4 (fc,q, ±l)m 5 (fc, q, ±l)m 6 (fc, q, ±1), 


dp 




= Tm^k, q, ±l)m 5 (k, q, ±l)m 6 (k, q, ±1), 
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L+4g 2 )/6 



Cf(*,q) = T m^(^q,±l)(e^- 

-p e -i(fe+2 ?1 +4 92 )/6^ 

Cf(fc,q) = T m 3 (A;,q,±l)m 6 (A;,q,±l)(e^ + ^ +4 '' 2 )/ 6 



=F e 



-i(fc+2qi+4g 2 )/6 



Cf } (/c,q) = =Fm 2 (A;,q, ±l)m 4 (fc,q, ±l)m 6 (fc, q, ±1), 



(182) 



and its dual. 
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